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The search for weakly interacting massive particles (WIMPs) by direct detection faces an en¬ 
croaching background due to coherent neutrino-nucleus scattering. As the sensitivity of these ex¬ 
periments improves, the question of how to best distinguish a dark matter signal from neutrinos will 
become increasingly important. A proposed method of overcoming this so-called “neutrino floor” 
is to utilize the directional signature that both neutrino and dark matter induced recoils possess. 

We show that directional experiments can indeed probe WIMP-nucleon cross-sections below the 
neutrino floor with little loss in sensitivity due to the neutrino background. In particular we And 
at low WIMP masses (around 6 GeV) the discovery limits for directional detectors penetrate be¬ 
low the non-directional limit by several orders of magnitude. For high WIMP masses (around 100 
GeV), the non-directional limit is overcome by a factor of a few. Furthermore we show that even for 
directional detectors which can only measure 1- or 2-dimensional projections of the 3-dimensional 
recoil track, the discovery potential is only reduced by a factor of 3 at most. We also demonstrate 
that while the experimental limitations of directional detectors, such as sense recognition and finite 
angular resolution, have a detrimental effect on the discovery limits, it is still possible to overcome 
the ultimate neutrino background faced by non-directional detectors. 

PACS numbers: 95.35.+d; 95.85.Pw 


I. INTRODUCTION 

Various cosmological observations indicate that ~ 30% 
of the energy density of the Universe is in the form of cold, 
non-baryonic, dark matter (CDM) [1]. Weakly Interact¬ 
ing Massive Particles (WIMPs) are a good CDM candi¬ 
date; they arise in extensions of the Standard Model of 
Particle Physics, such as Supersymmetry, and are nat¬ 
urally produced in the early Universe with the correct 
abundance (for reviews see e.g.. Refs. [H [3])- WIMPs 
from the dark matter halo of the Milky Way can be de¬ 
tected directly on Earth, via the keV-scale recoils pro¬ 
duced when they elastically scatter off nuclei [4]. 

Current direct detection experiments, including 
CDMS [3, LUX [5] and XenonlOO [7], are sensitive to 
spin-independent WIMP-nucleon cross-sections of order 
(Jx-n ~ 10“"^^ — 10“^® cm^. Signihcant increases in sensi¬ 
tivity are expected in the next few years as detector tar¬ 
get masses are increased to the ton-scale and beyond (see 
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e.g.. Ref. 0). As anticipated in early work on direct de¬ 
tection [3, these large detectors will also be able to de¬ 
tect coherent neutrino-nucleus scattering of astrophysical 
neutrinos [TOHT^ . Neutrinos are therefore the ultimate 
background for WIMP direct detection searches as they 
cannot be shielded against and produce recoils with sim¬ 
ilar rates and energy spectra [IIHIl]- 

For near-future direct detection experiments the most 
problematic types of neutrino are those produced in ®B 
decay in the Sun and in cosmic ray collisions in the 
Earth’s atmosphere. In a Xenon detector the recoil en¬ 
ergy spectrum and rate from ®B neutrinos very closely 
matches that of a WIMP with mass = 6 GeV and 
cross-section cr^_„ ~ 5 x cm^, while the spectrum 

from atmospheric neutrinos is similar to that of a WIMP 
with ~ 100 GeV and ~ 10“^® cm^ |T^. Gon- 
sequently the sensitivity of an experiment to WIMPs 
reaches a point of saturation where it becomes difficult 
to tell the difference between WIMP and neutrino in¬ 
duced recoils using their energies alone. So as the ex¬ 
posure or mass of an experiment increases the minimum 
discoverable cross-section rather than decreasing reaches 
a plateau. The value of cross-section at which this oc¬ 
curs depends on the systematic uncertainty in the neu- 



2 


trino flux and is commonly referred to as the “neutrino 
floor” [14]. It is worth noting however that for very 
large exposures, the neutrino floor is eventually mitigated 
when the small differences in the tails of the recoil en¬ 
ergy distributions of WIMPs and neutrinos start to dis¬ 
tinguish the two. However these required exposures are 
prohibitively large and well beyond the next generation of 
experiments so for practical purposes the neutrino back¬ 
ground does effectively present a “floor” to the discovery 
of dark matter and we will refer to it as such. 

If we wish to probe cross-sections below the neutrino 
floor then it is crucial to search for ways to distinguish the 
WIMP and neutrino signals, for instance via their differ¬ 
ent time and direction dependences. Grothaus et al. |15j 
have explored the sensitivity of directional detectors ca¬ 
pable of measuring the recoil directions in 3-dimensions 
using a hypothesis test which fits the WIMP and neu¬ 
trino event rates as a function of energy, time, and event 
angle with respect to the Earth-Sun direction. Davis 
[S] found that with very large exposures adding tim¬ 
ing information allows the neutrino floor to be evaded at 
low WIMP masses, due to the (small) annual modulation 
of both the WIMP and Solar neutrino signals. Ruppin 
et al. m examined how combining data from detectors 
composed of different target materials can probe cross- 
sections below the neutrino floor of a single experiment. 
They found that for spin-independent interactions the 
similarity in how the WIMP and neutrino signals scaled 
with respect to different target nuclei limited the advan¬ 
tage gained from multiple experiments. However for spin- 
dependent interactions the complementarity of multiple 
targets greatly improves the potential discovery limits. 

Directional detection experiments aim to reconstruct 
the nuclear recoil tracks in 3-dimensions. However this is 
experimentally challenging, and greater sensitivity might 
be achieved with a larger detector which measures the 1- 
or 2-dimensional projection of the recoil tracks. In this 
paper we extend the work of Ref. m by studying the 
effect of neutrino backgrounds on the sensitivity of ideal 
directional detectors with 1-, 2- and 3-dimensional read¬ 
out. We compare the potential discovery limits of these 
detectors with those of non-directional experiments that 
only measure the energy of the recoils, or only count the 
number of events above some threshold energy. We also 
extend our study beyond the ideal detector case and con¬ 
sider the effects of finite angular resolution and limited 
sense recognition. 

The paper is organized as follows. In Section |TT| we 
discuss the direction dependence of the cosmic neutrino 
fluxes. We then in Section lHll use these fluxes to calculate 
the directional neutrino event rates, and also review the 
calculation of the WIMP directional event rate. In Sec¬ 
tion [l^ we outline our analysis methodology for calculat¬ 
ing discovery limits. We present our results in Section [Vj 
Finally we conclude in Section |VI| with a summary and 
discussion of our results. 


II. NEUTRINO FLUXES 

In this section we review the Solar, atmospheric and 
diffuse supernovae background neutrino (DSNB) fluxes 
that dark matter detectors will be sensitive to. We ex¬ 
pand upon previous results in the literature by highlight¬ 
ing the angular dependence of these fluxes. 


A. Solar neutrinos 

Neutrinos produced from several reactions in the So¬ 
lar interior have been well-measured (for recent reviews 
see Refs. Most recently, the Borexino exper¬ 

iment has made the first spectral measurement of the 
pp component of the Solar neutrino flux [20] . The theo¬ 
retical systematic uncertainties on different components 
of the Solar neutrino flux range from 1% (pp flux) to 
14% (®B flux). For all measurements except the pp com¬ 
ponent, the theoretical uncertainties are as large as or 
larger than the measurement uncertainties. The theo¬ 
retical uncertainty arises largely from the uncertainty in 
the Solar metallicity, and in order to establish a self- 
consistent model of Solar neutrino fluxes one must as¬ 
sume a metallicity model. For the Solar neutrino flux we 
utilize the high metallicity Standard Solar Model (SSM) 
as defined in Ref. m-: low metallicity SSMs can predict 
lower flux normalizations by up to ~ 14% depending on 
the flux component. We consider here the high metallic¬ 
ity SSM because at present this model is more consistent 
with both the SNO neutral current measurement and he- 
liosiesmology data. In fact, future dark matter detection 
experiments will shed further light on the Solar metallic¬ 
ity issue [H]. Figure shows the fluxes and normaliza¬ 
tions of neutrinos produced from the different reactions. 
Due to their rather low energy (0(1 — 10) MeV) compared 
to the DSNB and atmospheric neutrinos. Solar neutrinos 
will mostly impact the discovery potential of future di¬ 
rect detection experiments in the low-mass region, below 
10 GeV (see Ref. [TH IT7] b 

In addition to the overall fluxes of Solar neutrinos, we 
are interested in their direction and time dependence. 
Due to the eccentricity of the Earth’s orbit, the Earth- 
Sun distance has an annual variation inducing a modula¬ 
tion in the Solar neutrino flux as seen by an Earth-based 
experiment such that, 

X — fl + 2ecos ( 

dFj^dH^ydt dEi, At \ T,y ) 

X d(q^ - q0(t)) , (1) 

where t is the time from January 1st, e = 0.016722 is 
the eccentricity of the Earth’s orbit, G = 3 days is the 
time at which the Earth-Sun distance is minimum (and 
hence the Solar neutrino flux is largest), = 1 year, At 
is the duration of the measurement, q,^ is a unit vector 
in the direction of interest and qo(t) is a unit vector in 
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FIG. 1: Neutrino energy spectra which are backgrounds to di¬ 
rect detection experiments: Solar, atmospheric, and the dif¬ 
fuse supernovae background. The Solar neutrino fluxes are 
normalized to the high metallicity SSM. The atmospheric neu¬ 
trinos are split into electron, antielectron, muon and antimuon 
neutrino components. The three DSNB spectra are labelled 
by their temperature in MeV, see Sec|II C| 


the inverse of the direction of the Sun^. As shown in 
Ref. [TS] , both the Solar neutrino and WIMP event rates 
have a ~ 5% annual modulation but they peak at times 
that are separated by about 5 months, and consequently 
timing information could help discriminate WIMPs from 
neutrinos. 


B. Atmospheric neutrinos 

At higher nuclear recoil energies, greater than approxi¬ 
mately 20 keV, the neutrino floor at high WIMP masses, 
i.e., above 100 GeV, will mostly be induced by low- 
energy atmospheric neutrinos (see [Bin]). These will 
limit the sensitivity of dark matter detectors without di¬ 
rectional sensitivity to spin independent cross-sections 
greater than approximately 10“^® cm^ [nmiiiiT]. 

The low energy flux of atmospheric neutrinos, less than 
approximately 100 MeV, is difficult to directly measure 
and theoretically predict [22]. At these energies, the un¬ 
certainty on the predicted atmospheric neutrino flux is 
approximately 20% [23]. Due to a cutoff in the rigidity 
of cosmic rays induced by the Earth’s geomagnetic field 
at low energies, the atmospheric neutrino flux is larger 
for detectors that are nearer to the poles P5] . 


^ We ignore the angular size of the Sun’s core on the sky which 
would give a tiny angular spread in the incoming neutrino direc¬ 
tions 


Over all energies, the atmospheric neutrino flux peaks 
near the horizon, at zenith angle cos 0 ~ 0. At high en¬ 
ergies, the flux is very nearly symmetric about cos d ~ 0, 
as at these energies the cosmic ray particles are more 
energetic than the rigidity cutoff. At low energies, the 
flux becomes asymmetric, as the flux of downward-going 
{cos 6 = 1) neutrinos is lower than the flux of upward- 
going neutrinos {cos 9 = —1). For the analysis in this 
paper, we consider the FLUKA results for the angular 
dependence of the atmospheric neutrino rate |24j . As we 
discuss below, we find that when this flux is convolved 
with the angular dependence of the coherent neutrino- 
nucleus cross-section, the angular dependence is washed 
out and the recoil spectrum depends only weakly on di¬ 
rection. There is also a seasonal variation in the neutrino 
flux based on the atmospheric temperature which induces 
an additional time modulation. However the exact time 
dependence of this effect at the latitude of our mock ex¬ 
periment is not known and is likely too small to have a 
large effect on the observed limits. Hence for this study 
we ignore both the angular and time dependence of the 
atmospheric neutrino flux and model it as isotropic and 
constant in time, 

1 

dEjydfl^dt An At dE^, 

C. Diffuse supernova neutrinos 

For WIMP masses between 10 and 30 GeV, the neu¬ 
trino floor is likely induced by the sub-dominant dif¬ 
fuse supernova neutrino background (DSNB), from all 
supernova explosions in the history of the Universe. The 
DSNB flux is a convolution of the core-collapse supernova 
rate as a function of redshift with the neutrino spectrum 
per supernova; for a recent review of the predicted DSNB 
flux see Beacom j^Sj. The DSNB spectra have a similar 
form to a Fermi-Dirac spectrum with temperatures in 
the range 3-8 MeV. We use the following temperatures 
for each neutrino flavour: = 3 MeV, = 5 MeV 

and = 8 MeV, where Vx represents the four remaining 
neutrino flavours. Motivated by theoretical estimates we 
take a systematic uncertainty on the DSNB flux of 50%. 
The DSNB is believed to be isotropic and constant over 
time, therefore its angular dependence can be expressed, 
as with the atmospheric neutrinos, using Eq. ([^. 

III. NEUTRINO AND DARK MATTER RATE 
CALCULATIONS 

A. Coherent neutrino-nucleus elastic scattering 

We only consider the neutrino background from coher¬ 
ent neutrino-nucleus elastic scattering (GNS) as it pro¬ 
duces nuclear recoils in the keV energy scale which cannot 
be distinguished from a WIMP interaction. We neglect 
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1 / type 

(MeV) 

(keV) 

u flux 

(cm-®s-®) 

PP 

0.42341 

2.94 X 10"^ 

(5.98 ±0.006) X 10“ 

■'Be 

0.861 

0.0122 

(5.00 ±0.07) X 10® 

pep 

1.440 

0.0340 

(1.44 ±0.012) X 10® 

ISq 

1.732 

0.04917 

(2.23 ±0.15) X 10® 

«B 

16.360 

4.494 

(5.58 ±0.14) X 10® 

hep 

18.784 

5.7817 

(8.04 ± 1.30) X 10® 

DSNB 

91.201 

136.1 

85.5 ±42.7 

Atm. 

981.748 

15.55 X 10® 

10.5 ±2.1 


TABLE I: Dominant neutrino fluxes with corresponding un¬ 
certainties. For the Solar neutrino flux, we utilize the high 
metallicity SSM, as described in the text. The maximum neu¬ 
trino energy, and maximum recoil energy on a Xenon 

target, are also shown. 


neutrino-electron elastic scattering, mostly induced by pp 
neutrinos, as it has been shown to only marginally affect 
the discovery potential of experiments with limited nu¬ 
clear/electronic recoil discrimination for WIMP masses 
above 100 GeV [Ti] . 

Freedman [26] has shown that neutrino-nucleon elastic 
scattering, which is well explained by the standard model 
but has yet to be observed, leads to a coherence effect at 
low momentum transfer that approximately scales with 
the atomic number of the target nucleus. A, squared. 
At higher recoil energies, generally above a few tens of 
keV, the loss of coherence is described by the nuclear 
form factor F{Er), for which we use the standard Helm 
form m- The differential cross-section as a function of 
the nuclear recoil energy {Er) and neutrino energy {E^) 
is given by 


setting /3 = 0 in Eq. Q, 

^max ^ _ 2El 

^ [Ei,+mNY mN + ‘2.E^ 


( 6 ) 


The maximum recoil energies produced by the different 
types of neutrino for a Xenon target are shown in Tablejl] 
The directional event rate per unit mass and time, as a 
function of the recoil energy, direction and time, is given 
by the convolution of the double differential CNS cross- 
section and the neutrino directional flux. 


d'^i? 


dErdHr-dt 




dV 

dErdflr 


d'^$ 


(7) 

where A™'" = •^miqEj./2 is the minimum neutrino en¬ 
ergy required to generate a nuclear recoil with energy E^ 
and .yV is the number of target nuclei per unit mass. 

The directional event rate for Solar neutrinos is found 
by substituting Eqs. 0 , 0 and 0 into Eq. 0 and 
integrating over the neutrino direction 


d?R 


dErd^r<it 
f dcr dd) 

J dEr dE^ 


27r 


S 


At 

1 ± 2e cos 

q© - 

El, ± Tow 

El, 


/ 27r(t - t^) 

V 


Er 


2m N 


dE,. (8) 


The delta function can then be rewritten as 


j . Eu + ruN 

I Qr • q©-- 



1 

Z?mirL 


where we have defined x = — IjEu and, 


(5 



(9) 


1. ^ qr • q© _1_ 

£ A™" TOW ■ 


( 10 ) 


where Qw = M — — 4sin^ &w)Z is the weak nuclear 

hypercharge of a nucleus with M neutrons and Z pro¬ 
tons, Gf is the Fermi coupling constant, 6^ is the weak 
mixing angle and tow is the target nucleus mass. The 
directional and energy double differential cross-section 
can be written by noting that the scattering has az¬ 
imuthal symmetry about the incoming neutrino direc¬ 
tion so dHj^ = 27r d cos jd and imposing the kinematical 
expression for the scattering angle, /3, between the neu¬ 
trino direction, q,^, and the recoil direction, q^, 


cos [i = cir- hu = 


E^ + / E, 


E. 


2tow 


(4) 


with fd in the range (0,7r/2), using a delta function. 


d (T 


dErdHr dEj. 27r 


dcr 1 - / ^ E^ + mN I Er 

— 0 cos p — 


E„ 


2m N 


(5) 

The maximum recoil energy, can be obtained by 


Finally, by substituting Eq. 0 into Eq. 0, integrat¬ 
ing over X and converting back to E^, we obtain an ana¬ 
lytic expression for the directional event rate from Solar 
neutrinos, 


d^R 

dErdHj.dt 


27r At 


1 -I- 2e cos 


/ 27r(t - U) 
\ 


X 


£'^ dcr 

l^dlA 


{Er,£) 


d$ 

di^ 


e 


( 11 ) 


for cos“^(qr. • q©) < tt/2 and 0 otherwise. 

In the case of the atmospheric and diffuse supernova 
neutrinos, as we have assumed their fluxes to be isotropic 
and constant over time (see Sec.|0, the directional event 
rate is simply given by substituting Eqs. 0 and 0 into 
Eq. 0 and integrating over the neutrino direction 
leading to 


d^R _ wF f ^ —dE 

dErdflj-dt 47rAt dEr 'dE^, ' 


( 12 ) 
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FIG. 2: Neutrino and = 6GeV WIMP nuclear recoil rates for a Xenon target as a function of recoil energy, Er, (left) and 
cosine of the angle between the Solar vector and recoil vector, cos^sun, (right) obtained by integrating the differential recoil 
spectrum over angle and energy respectively. The atmospheric and DSNB neutrino fluxes are taken to be isotropic. In order 
to show all types of neutrino the threshold energy has been set to 1 eV here. 


The neutrino event rates as a function of energy and 
angle between Solar and recoil directions, cos^sun = 
—Qr ■ q0, obtained by integrating Eqs. ( [Tl] ) and (12) 
over direction and energy respectively are shown in Fig.[^ 
Also shown is the recoil rate for a 6 GeV WIMP, showing 
the similarity between this spectrum and the spectrum 
of ®B neutrino recoils. The isotropic DSNB and atmo¬ 
spheric recoil rates are flat whereas the event rates of 
Solar neutrinos are highly anisotropic. The curves corre¬ 
sponding to the mono-energetic neutrinos (^Be and pep) 
have a sharp cutoff in their directionality due to the finite 
energy threshold. From Fig. [^one can already anticipate 
that the degeneracy between solar neutrino and WIMP 
events from an energy-only analysis will be almost com¬ 
pletely removed with the addition of directional informa¬ 
tion. 


B. Dark matter 

Like most spiral galaxies, the Milky Way is believed to 
be immersed in a halo of dark matter which outweighs 
the luminous component by at least an order of mag¬ 
nitude P51l30| . The velocity distribution of dark mat¬ 
ter in the halo is traditionally modelled (cf. Ref. E3) 
as an isotropic Maxwell-Boltzmann distribution, corre¬ 
sponding to an isotropic halo with a I/r^ density profile 
and a flat rotation curve. Simulated halos have veloc¬ 
ity distributions which deviate systematically from the 
Maxwell-Boltzmann distribution [3TV137) . However, these 
deviations do not have a large effect on the discovery 
potential of directional detection experiments [3S1 155] . 


Therefore, and to allow comparison with previous work, 
we assume a truncated Maxwell-Boltzmann distribution 
throughout this work which in the Galactic rest frame 
has the form, 


/gai(v) 


Aresc(27I-<T2)3/2 

1 

1 

0 


if |v| < Vesc , 
if |v| > Vesc , 


(13) 

where ay is the WIMP velocity dispersion which is related 
to the local circular speed, vq, by ay = vojV^, Vesc is the 
escape speed, and A^esc is a normalization constant, 


Nesc = erf 


'^esc \ 

V^ayJ 



(14) 


We use the conventional values for the circular and escape 
speeds: vq = 220kms“^ [40| and Vesc = 544kms“^ [4T] . 

The velocity distribution of WIMPs in the rest frame 
of the laboratory is obtained through a Galilean trans¬ 
formation of the Galactic frame distribution, /gai, by the 
laboratory velocity Viab (discussed below). 


/lab(v) = /gal(v -b Viab) . (15) 


The directional event rate as a function of both recoil 
energy, direction in the lab frame and time, assuming 
spin-independent interactions with identical couplings to 
protons and neutrons, is given by [l^l 155] 


d^R 




dEydflydi 47rm^/:i^„At 




n Qri i): 

(16) 
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respect to the local standard of rest, vsoiar, the Earth’s 
revolution around the sun, VEarthRev, and the Earth’s ro¬ 
tation, VEarthRot- The dominant contribution to the lab 
velocity is the sum vcaiRot + vsoiar while VEarthRev and 
VEarthRot are, respectively, responsible for the annual [5] 
and diurnal m modulation effects. Even though their 
contributions to Viab are small, we take them into account 
in order to accurately study the additional discrimination 
power that can be brought by the annual and diurnal 
modulation effects. A detailed review of the Galactic- 
to-lab frame velocity and coordinate transformations are 
given in Ref. m- In the following, we will consider a 
detector with the x, y, and z axes pointing toward the 
North, the West and the Zenith directions respectively. 
For a detector located at a latitude 90°, i.e., at the North 
Pole, z is aligned with the spin axis of the Earth. Finally, 
in the detector frame, the direction of a recoil is given by 
the angles 9 and (j) defined such that. 


FIG. 3: The position in the sky, in terms of Altitude and Az¬ 
imuth, of the Sun (red) and the inverse of the position of the 
constellation Cygnus (corresponding to the direction -viab) 
(blue) as observed from the Modane underground laboratory 
(latitude 45.2°). The points show the position at measure¬ 
ments made every hour from the 1st January 2015 0:00 until 
31st December 2015 23:00. The Solar position traces out 24 
analemmas, corresponding to the Sun’s position at each hour 
of the day over the course of a year. The dashed horizontal 
line is the horizon. As demonstrated here, the Sun’s position 
does not coincide with that of Cygnus at any time. 


where is the WIMP mass, the WIMP-nucleon 
reduced mass, po = 0.3 GeV cm“^ the local dark mat¬ 
ter density, A the mass number of the target, cFx-n the 
WIMP-nucleon cross-section and Vmin = \/2 to)vE/2p^„ 
is the minimum WIMP speed required to produce a 
nuclear recoil of energy E^. Finally, /lab (f min j Qr; t) 
is the three-dimensional Radon transform of the lab- 
frame WIMP velocity distribution /lab(v), which for the 
Maxwell-Boltzmann distribution has the form [43] . 


/lab('Cmin7 Qr! t) 


X 



^■380(2^2) 1/2 

Imin + q ■ Viab I ^ 
2cr2 


— exp 



As can be seen from this equation, if Umin < iiiab the 
directional rate is maximum for q • vjab = —Vmin (i-S., 
in a ring [44] ). while otherwise the rate is maximum for 
q = —Viab (he., a dipole distribution [H]). The strong 
correlation between the recoil directions and the labora¬ 
tory motion in the Galactic frame allows the unambigu¬ 
ous authentication of a WIMP signal [H] [46| . 

The lab velocity, viab, is given by the sum of the ro¬ 
tation of the Solar System around the Galactic center 
VGaiRot, the peculiar velocity of the Solar System with 


q = sin 6 cos (j)yi + sin 9 sin (j>y + cos 9 z . (18) 

The laboratory velocity viab happens to point towards 
the constellation of Gygnus. Figure [^ shows the position 
in the sky of the Sun and the inverse of the position of 
Cygnus (i.e., —viab), as observed from the Modane un¬ 
derground laboratory (latitude 45.2°). The points show 
the position at observations made every hour from the 
1st January 2015 0:00 until 31st December 2015 23:00. 
The Solar position traces out 24 analemmas, correspond¬ 
ing to the Sun’s position at each hour of the day over the 
course of a year. As we see here, the Sun’s position does 
not coincide with that of Cygnus at any time suggest¬ 
ing that a directional experiment should in principle be 
able to disentangle the WIMP from the Solar neutrino 
contributions in the observed data. As a matter of fact, 
the angular separation between the peak WIMP direction 
and the peak neutrino direction undergoes a sinusoidal 
modulation over the course of the year that varies from 
60° in February to 120° in September. 


C. Resulting signals 

Figure [^ shows Mollweide projections of the labora¬ 
tory frame angular differential event rate from a 6 GeV 
WIMP plus ®B Solar neutrinos, at the times when the 
separation between the directions of the Sun and Cygnus 
are smallest (60°) and largest (120°). This Figure clearly 
shows that, even at the time of smallest separation, the 
WIMP and neutrino recoil distributions can be easily dis¬ 
tinguished as long as the angular resolution is better than 
a few tens of degrees. Although this Figure only shows 
the rates for ®B neutrino induced recoils, the angular 
distributions for other Solar neutrinos are very similar as 
neutrinos can only induce a recoil with an angle in the 
range (0,7r/2) from their incident direction. Additionally 
the angular dependence of the recoil spectra is correlated 
with energy as can be seen going from the left to the right 
hand panels. For both the WIMP and neutrino recoils 
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Time: 26th Feb. 2015 06:00 Target: Xe WIMP: = & GeV, = 4.9 x 10^^® cm^ 

0 - 1.6667 keV 1.6667 - 3.3333 keV 3.3333 - 5 keV 



0 100 200 300 400 0 5 10 15 20 0.01 0.02 0.03 0.04 

Time: 6th Sep. 2015 06:00 


0 - 1.6667 keV 1.6667 - 3.3333 keV 3.3333 - 5 keV 



0 100 200 300 0 5 10 0 0.01 0.02 0.03 0.04 

di^bin/'^^r year^^ sr^^] 


FIG. 4: Mollweide projections of the WIMP plus ®B neutrino angular differential event rate integrated within (from left 
to right) three equally sized energy bins spanning the range Er = Q to 5keV, for a WIMP with mass = 6 GeV and 
(Jx-n = 4.9 X cm^ and a Xe target. The top row shows the signal on February 26th, when the separation between the 

directions of the Sun and Gygnus is smallest (~ 60°), and the bottom row on September 6th, when the separation is largest 
(~ 120°). The WIMP contribution is to the left of the neutrino contribution on the top row and to the right on the bottom 
row. The Mollweide projections are of the event rate in the laboratory co-ordinate system with the horizon aligned horizontally 
and the zenith and nadir at the top and bottom of the projection respectively. 


the angular spread decreases with increasing energy i.e., 
the highest energy recoils have the smallest angle between 
the incoming particle direction and the recoil direction. 


In addition to the standard case of a detector with full 
3-dimensional sensitivity, we will also assess the discov¬ 
ery potential of a detector which only has sensitivity to 
1-dimensional and 2-dimensional projections of the 3-d 
recoil track. Using Eq. (18) we define the 2-d readout to 
be the projection of the recoil track onto the x-y plane 
such that only the angle 4> is measured, and the 1-d read¬ 
out to be the projection on to the x-axis such that only 
the angle 6 is measured. 


Figure shows the daily evolution of the 1-d, cos0, 
and 2-d, (/>, recoil angle distributions at a single energy 
(0.5 keV) from ®B neutrinos and a WIMP with mass 
TOjj. = 6 GeV. The (f> distributions from ®B neutrinos 
have two peaks, because at a fixed recoil energy the neu¬ 
trino energy spectrum produces recoils in a ring around 
the incident direction. In the WIMP case, however, the 
distribution of recoils is peaked in a single direction, to¬ 
wards —viab. The 2-d and 1-d distributions for both at¬ 


mospheric and DSNB neutrinos are flat, and therefore we 
do not show them for clarity. The WIMP and neutrino 
distributions are significantly different, not only in their 
shape at a single time but also how they evolve over the 
course of a day. This suggests that a detector with only 
1-d or 2-d readout should still be able to discriminate 
WIMP and neutrino induced recoils. 


IV. ANALYSIS METHODOLOGY 

In this section we introduce the analysis methodology 
we use to assess the discovery potential of each readout 
strategy for future low-threshold, ton-scale experiments. 
Discovery limits were first introduced in Ref. |39] and are 
defined such that if the true WIMP model lies above this 
limit then a given experiment has a 90% probability to 
achieve at least a 3a WIMP detection. To derive these 
limits, it is necessary to compute the detection signif¬ 
icance associated with different WIMP parameters, for 
each detector configuration. This can be done using the 

















FIG. 5: The daily evolution, at three hourly intervals, of the angular distributions of 0.5 keV Xe recoils from ®B neutrinos (left 
column) and a WIMP with mass = 6GeV (right). The distributions are normalized to unity in each case and displayed 
with arbitrary units. The top row shows the distribution of cos 8 measured by a detector with 1-d readout and the bottom row 
the angle ij) for 2-d readout. See the text for the dehnitions of 9 and 4>- The date chosen was September 6th 2015, the date of 
maximum separation between the WIMP and neutrino distributions. 


standard profile likelihood ratio test statistic [48] where 
the likelihood function at a hxed WIMP mass is defined 
as 




= -(Mx+Ej = i tJ'i) 


Nl 


N 

n 


Mx/x (q* > ) 


i=i 




(19) 


k=l 


where the sums j and k are over the n,^ neutrino back¬ 
grounds, and N are, respectively, the expected 

number of WIMP and neutrino events, and the total 
number of observed events, and fl are the normal¬ 
ized, time and momentum dependent, event rates for the 
WIMP and neutrinos, ti is the time at which the event 
occurred and corresponds to the set of observables for 
each event, which depends on the readout considered. 
For 3-d readout = {£^,0,(1)}. Finally, are 

the individual likelihood functions associated with the 


flux <I>fe of each neutrino component. These individual 
likelihood functions are each parametrized as gaussian 
distributions with a standard deviation given by the rel¬ 
ative uncertainty in the neutrino flux normalization as 
discussed in Sec. |n] (see Table |I]). 

The profile likelihood ratio corresponds to a hypothesis 
test between the null hypothesis Hq (background only) 
and the alternative hypothesis Hi which includes both 
background and signal, incorporating systematic uncer¬ 
tainties in this case the normalization of the neutrino 
fluxes. As we are interested in the WIMP discovery po¬ 
tential of future experiments, we test the background 
only hypothesis, Hq, on simulated data and try to re¬ 
ject it using the following likelihood ratio. 


A(0) 


= 0 ,^) 


( 20 ) 


where $ and denotes the values of $ and that 

maximize the unconditional .jSf and # denotes the values 
of $ that maximize .if under the condition (T^_n = 0, 
i.e., we are profiling over the parameters in $ which are 
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considered to be nuisance parameters. As discussed in 
Ref. |48], the test statistic qo is then defined as, 


qo = 


—21nA(0) ay.-n > 0 , 
0 ^x-n ^ 0 - 


( 21 ) 


and that the 9 and (j) angles are defined as in Eq. (181. 
Following previous work dedicated to the comparison of 
different readout strategies in the context of an arbitrary 
background [SI], we consider 6 detector readout strate¬ 
gies: 


A large value for this statistic implies that the alterna¬ 
tive hypothesis gives a better fit to the data, i.e., that it 
contains a WIMP signal. The p-value, po, of a particu¬ 
lar experiment is the probability of finding a value of qo 
larger than or equal to the observed value, if the 

null (background only) hypothesis is correct: 

POO 

Po= f{qo\Ho)dqo, (22) 

where f{qo\Ho) is the probability distribution function 
of qo under the background only hypothesis. Following 
Wilk’s theorem, qo asymptotically follows a distribu¬ 
tion with one degree of freedom (see Ref. |48| for a more 
detailed discussion) and therefore the significance Z in 
units of standard deviation is simply given by Z = \/q^- 
The discovery limit for a particular input WIMP mass 
can then be found by finding the minimum cross-section 
for which 90% of the simulated experiments have Z > 3. 
In the following results we have used distributions from 
5000 Monte-Carlo experiments. 

This analysis methodology was first introduced by the 
XENONIO collaboration |7] and many experiments are 
now using similar likelihood approaches e.g., LUX [Bj, 
CDMS-II [49||50], and CoGeNT [51]. This has become 
possible thanks to the construction of accurate back¬ 
ground models derived from reliable simulations, as well 
as data-driven analysis techniques based on calibration 
data. The advantage of using likelihood analyses is that 
they can not only determine whether or not a dark mat¬ 
ter interpretation to the data is preferred and a WIMP 
signal detected, but they can also measure or constrain 
the WIMP parameters themselves. Furthermore a like¬ 
lihood analysis will maximize the sensitivity to the dark 
matter signal and obtain the best possible limits for a 
given experiment. 

V. RESULTS 
A. Detector configurations 

As the goal of this paper is to give a detailed overview 
of how the neutrino background will affect future di¬ 
rect detection experiments, we first give a brief descrip¬ 
tion of the different readout strategies that we consider. 
Here, and throughout, we consider a Xe-based exper¬ 
iment located in Modane with latitude and longitude 
(45.2°, 6.67°), taking data over a duration At = 1 year. 
As mentioned above, we assume that the reference frame 
of the detector is such that x, y, and z are, respectively, 
pointing toward the North, West, and Zenith directions 


• 3-d directional readout, {Er,9,(j),t} 

• 2-d directional readout, {Er,(^,t} 

• 1-d directional readout, {Er,9,t} 

• No directional information {Er,t} 

• Event time only {t} 

• Number of events only (i.e., a counting experiment) 

The last two strategies correspond to detectors that 
can only measure the total number of events above some 
threshold. This is the case for bubble chamber exper¬ 
iments [53] that adjust their operating pressure to nu¬ 
cleate a single bubble from a nuclear recoil. The energy 
and time, {E^, t}, strategy corresponds to the majority of 
current and ongoing direct detection experiments where 
the kinetic energy of the recoiling nuclei is obtained from 
measurements of the heat, ionization and/or scintillation 
energies deposited in the detector (see Ref. |S3] for a re¬ 
cent review). 

The 3-d directional readout, {Er,9,(p,t}, corresponds 
to the ultimate detector that measures and exploits all 
the information available in the WIMP recoils. Current 
directional experiments are using low-pressure gaseous 
Time Projection Chambers (TPCs) in order to obtain 
tracks from 0(10) keV nuclear recoils that are a few 
mm long (see Ref. (54] and references therein). For 3- 
d sensitive directional detectors, the track is measured 
by sampling over time the 2-dimensional projection of 
the ionization-induced electron cloud on a pixelized an¬ 
ode. For a 2-d readout, {Er,(j),t}, which is usually 
based on CCD technology, the anode is not time-sampled 
and therefore only a 2-dimensional projection of the 
drifted electron cloud can be measured. A 1-d read¬ 
out, {Er,9,t}, only measures the projection of the re¬ 
coil track along the drift direction. This could be done 
for example in a dual-phase Liquid Xe TPC experiment 
by looking at the ratio of the ionization and scintillation 
energies thanks to columnar recombination of the drift¬ 
ing electrons/ions |SS|. However, this effect has yet to 
be confirmed and ongoing measurements by the SCENE 
collaboration have only found mild evidence for columnar 
recombination from keV-scale nuclear recoils [56] . 

B. Comparing readout strategies 

In this section we consider idealised detectors with per¬ 
fect efficiency and angular resolution and no backgrounds 
apart for neutrinos, hence the discovery limits we obtain 
represent the best-case scenario for a particular read¬ 
out strategy. We will study experimental limitations in 
Sec. IV Cl and IV PI 
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FIG. 6: The dependence of the discovery limit for the spin independent WIMP-nucleon cross-section, on the mass of a 

Xe detector operated for lyear using (from top to bottom) the number of events only (pink line), time information (brown 
dotted), energy & time (orange), energy & time plus 1-d (red), 2-d (blue) and 3-d (green) directionality. The left (right) plot 
is for = 6 (100) GeV and an energy threshold Sth = 0.1 (5) keV and the bottom axis shows the number of ®B (atmospheric) 
neutrinos expected. Note the different scales of the left and right hand plots. 


Figure 1^ shows the evolution of the spin-independent 
discovery limit with detector mass, M, for a Xe detector 
taking data for one year using each of the different read¬ 
out strategies. We consider two example WIMP masses 
and detector: a light WIMP & low threshold detector 
{m^ = 6 GeV and Eth =0.1 keV) and a 100 GeV 
WIMP & high threshold detector {Eth = 5 keV). For 
these two WIMP masses the recoil energy spectra closely 
matches those of ®B and atmospheric neutrinos respec¬ 
tively. As found previously, when the expected number 
of neutrino background events is negligible, the discovery 
limits improve rapidly with detector mass as a function of 
1/M [13] and the difference between the readout strate¬ 
gies is very small (cf. Ref. [52]). As the detector mass is 
increased and the experiment begins to have an apprecia¬ 
ble neutrino background a Poisson background subtrac¬ 
tion regime is entered and the discovery limit evolves as 
I/'/m. When the expected number of neutrino events 
reaches 10 — 10^ the counting only, time only, and en¬ 
ergy & time limits plateau at a value controlled by the 
systematic uncertainty on the dominant neutrino compo¬ 
nent according to El: 


CTDL OC 



(23) 


where udl is the discovery limit, ^ is the uncertainty on 
the relevant neutrino flux (Table|T| and /j,^ is the expected 
number of neutrinos. In this regime the experiment can¬ 
not tell the difference between WIMP and neutrino in¬ 
duced recoils as there are not enough events to probe the 
different time dependences, or the differences in the tails 


of the energy distributions. This saturation of the WIMP 
sensitivity, which spans over two orders of magnitude in 
exposure, is what is commonly referred to as the neutrino 
floor. 

The limits with directional readout however continue 
to decrease as the incorporation of directional informa¬ 
tion allows the distributions of WIMP and neutrino in¬ 
duced recoils to be distinguished. For the 100 GeV 
WIMP case, the limits from 2-d and 3-d readout are 
a factor of ~ 1.2 and 1.6 better than those from 1-d 
readout whereas for the 6 GeV WIMP case they are fac¬ 
tors of ^ 1.2 and 3 times better. The discovery limit 
with directionality continues to decrease as 1/M for the 
6 GeV WIMP as the directional and time-dependent dis¬ 
tributions of the WIMP and Solar neutrino induced re¬ 
coils have only very small overlap such that the back¬ 
ground has very little effect on the discovery capabilities 
of the experiment. However for the 100 GeV WIMP the 
dominant background from isotropic atmospheric neu¬ 
trinos significantly overlaps with the WIMP distribu¬ 
tion so, although the experiment is able to distinguish 
the WIMP signal, the sensitivity is still compromised by 
the background. Therefore, in this case the discovery 
limit, beyond the saturation regime, evolves according 
to a standard Poisson background subtraction mode as 
1/^/M. Our results are generally consistent with those of 
Grothaus et al. ng, in that we agree that directionality 
is the most promising strategy to go beyond the neu¬ 
trino floor. However, there are quantitative differences. 
We find that directionality allows greater improvements 
in sensitivity than found in Ref. m- This is largely 
due to the difference in energy thresholds used (2 keV 
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FIG. 7: The discovery limit as a function of WIMP mass 
using (from top to bottom) the number of events only (pink 
solid line), time information (brown dotted), energy & time 
(orange), energy & time plus 1-d (red), 2-d (blue) and 3-d 
(green) directionality. The upper (lower) set of lines are for 
the two detector set-ups described in Table|n] Detector A (B) 
with a target mass M — 0.1 (10^) ton and an energy threshold 
Sth = 0.1 (5) keV. The black curve and shaded region shows 
the neutrino floor from Ref. [la¬ 


in Ref. [TS] compared to 0.1 keV in this work) which 
drastically changes the ratio of WIMP to neutrino event 
numbers for low mass WIMPs. We also believe that some 
differences are a result of the fact that we have used the 
full directional information rather than just the 

reduced angles, and the two analyses also use different 
statistical techniques. We have hitherto considered an 
ideal detector, however finite angular resolution at the 
level considered in Ref. m does not significantly change 
our conclusions (see Sec. VD). 

For very large detector masses (M > 10 ton for Eth = 
0.1 keV and M > 10"^ ton for i?th = 5 keV) which have 
accumulated more than ~ 10“^ neutrino events, the evolu¬ 
tion of the Time only and Energy -|- Time discovery lim¬ 
its return to the Poisson background subtraction regime 
once more. With a very large number of events the time 
information allows discrimination between WIMP and 
neutrino induced recoils (cf. Ref. |16|). However time 
information is more useful for discriminating Solar neu¬ 
trinos from light WIMPs than for discriminating atmo¬ 
spheric neutrinos from heavier WIMPs. This is because 
the WIMP and Solar neutrino rates are both annually 
modulated, and also the amplitude of the annual modu¬ 
lation is larger for light WIMPs. For energy information 
only, with very large numbers of events the slight dif¬ 
ference in the tails of the ®B neutrino and WIMP recoil 
energy distributions allows them to be discriminated [HI- 

Having studied the evolution of the discovery limit as a 
function of detector mass for two specific WIMP masses. 


we now consider two fixed example detector set-ups out¬ 
lined in Table El a low mass & low threshold detector 
(M = 0.1 ton and Eth = 0.1 keV respectively) and a 
high mass & high threshold detector (lO'^ ton and 5 keV). 
Again, for simplicity and to probe the full annual modula¬ 
tion signal, we assume that data is accumulated over one 
year. These detector masses and thresholds are chosen so 
that a non-directional detector with the same mass and 
threshold would be in the saturation regime that results 
in the neutrino floor, as seen in Fig. 

Figure [7] shows the discovery limit as a function of 
WIMP mass for the two detector set-ups for each read¬ 
out strategy. Also shown as the shaded region is the neu¬ 
trino floor from Ref. m which is the combination of two 
limits obtained by a Xenon detector. For light WIMPs 
{m^ < 10 GeV) the limit comes from a 3 eV thresh¬ 
old detector with an exposure of 0.19 ton years, while 
for heavier WIMPs (m^ > 10 GeV) a detector with a 4 
keV threshold and an exposure of 9.3 x 10^ ton years was 
used. The two detector configurations roughly match our 
two detector setups A and B in Table |Hj As described in 
Refs. [131 [T7], the low-mass part of the neutrino floor 
comes from solar neutrinos (which have low-energies but 
high fluxes) with the shoulder at = 6 GeV arising due 
to ®B neutrinos. The high-mass part, above ^ 10 GeV, is 
due to DSNB neutrinos and atmospheric neutrinos which 
have higher energies but much lower fluxes. 

For the low threshold detector, the directional discov¬ 
ery limits clearly cut through the low-mass neutrino floor 
and for the 3-d readout there is actually almost no reduc¬ 
tion in sensitivity due to the neutrino background. The 
1-d and 2-d readouts do suffer a small reduction in sensi¬ 
tivity, but evidently the distributions are different enough 
that it is still possible to probe cross-sections below the 
limit set by non-directional experiments. For the high 
threshold detector the improvement in the discovery lim¬ 
its, with respect to the high-mass neutrino floor, from di¬ 
rectionality is smaller. However it does still help discrimi¬ 
nate the isotropic atmospheric neutrino background from 
WIMP induced recoils, in particular for WIMP masses 
around 100 GeV where the energy spectra from WIMPs 
and atmospheric neutrinos are most similar. 

In summary, we found that directionality is a pow¬ 
erful tool for disentangling neutrino backgrounds from 
a putative WIMP signal. The gain from directionality 
is particularly impressive for low mass WIMPs thanks 
to the large separation between the solar neutrino and 
WIMP incoming directions, see Sec. |HIC[ Interestingly, 
we found that this result still holds even if only the 2-d or 
1-d projection of the recoil tracks can be measured. The 
gain from directionality in the high-mass region is more 
moderate, however, due to the large overlap between the 
WIMP and the isotropic DSNB and atmospheric neu¬ 
trino distributions. Even in this case, we found that I-d 
and 2-d readouts still outperform non-directional experi¬ 
ments. This highlights that it is worthwhile to construct 
directional detectors, even without full 3-d readout. Sim¬ 
ilar conclusions were reached in Ref. [55] in the context 
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Detector 

Target 

M (ton) 

Eth (keV) 

At 

(Lat., Long.) 

1 / backgrounds 

A 

Xe 

0.1 

0.1 

1 yr 

(45.2°, 6.67°) 

®B, hep, Atm., DSNB 

B 

Xe 

10^ 

5 

1 yr 

(45.2°, 6.67°) 

Atm., DSNB, hep 

[IT] (low) 

Xe 

0.19 (ton-year) 

0.003 

- 

■^Be, 

®B, pep, ®®0, ®®N, ®^F, hep. Atm., DSNB 

m (high) 

Xe 

9.3 X 10® (ton-year) 

4 

- 

- 

Atm., hep, DSNB 


TABLE II: The properties of the two detector set-ups we consider: M is the target mass, i?th the energy threshold and At the 
length of time, from Jan. 1st 2015, over which data is taken. In both cases the detector is located at Modane. For reference we 
also show the properties of the detector set-ups used in Ref. m to generate the neutrino floor shown in Fig. and subsequent 
plots. Also listed are the neutrino backgrounds present in each experiment, in order of number of expected events from highest 
to lowest. 



FIG. 8: The discovery limit as a function of the energy 

threshold for head-tail discrimination for 1-d (red), 2-d (blue) 
and 3-d (green) directional readout (with energy and time in¬ 
formation in all three cases). The dashed black line shows the 
discovery limit with energy information only. The left (right) 
hand set of curves is for = 6 (100) GeV and detector set-up 
A (B), with low (high) mass and threshold. 



FIG. 9: The discovery limit as a function of WIMP mass 
for detectors with full sense recognition (solid lines) and no 
sense recognition (dashed) and 1-d (red), 2-d (blue) and 3-d 
(green) directional readout. The upper (lower) set of lines are 
for detector set-up A (B), with low (high) mass and threshold. 
The dashed black lines show our discovery limit with energy 
information only and the black curve and shaded region shows 
the neutrino floor from Ref. m- 


of a discriminating WIMPs from an arbitrary isotropic 
background. 


C. Sense recognition 

As well as the 3-dimensional reconstruction of the re¬ 
coil track, a key experimental concern is head-tail recog¬ 
nition i.e., the ability to measure the sense (-fq or — q) 
of the nuclear recoils. Determining the sense of a recoil 
is expected to be possible in a gas TPC by measuring 
the asymmetry in the charge collected along the recoil 
track as well as asymmetry in the shape of the track 
itself m- Whilst there has been important progress 
with the experimental effort, sense recognition remains 
one of the outstanding challenges for current and fu¬ 
ture detectors 133] . Whether or not this information 
can be retrieved from a detector plays a key role in its 


ability to discriminate WIMP induced recoils from back¬ 
grounds [38l [39l [60] . This is because in the absence of 
head-tail discrimination the angular recoil rates in the 
forward and background directions are added together 
and the anisotropy of the WIMP induced recoils is effec¬ 
tively decreased. 

Figure]^ shows how the discovery limit depends on the 
energy threshold for head-tail discrimination, for = 6 
and 100 GeV and our two example detector set-ups (Ta¬ 
ble [IT]) with 3-d, 2-d and 1-d readouts. For simplicity, we 
assume that above (below) the head-tail energy thresh¬ 
old there is perfect (no) head-tail discrimination. For 
the light WIMP and the low mass and threshold detec¬ 
tor, the discovery limits are weakened as the head-tail 
energy threshold is increased from 0.1 keV to ~ 1 — 2 keV 
before flattening off to a factor between ~ 1.5 (1-d) and 
~ 10 (3-d) below the energy only limit. For lower dimen- 
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FIG. 10: The dependence of the discovery limit for the spin independent WIMP-nucleon cross-section, on detector mass 

for detectors with (solid lines) and without (dashed lines) sense recognition for each of the three directional readout strategies 
3-d (green), 2-d (blue) and 1-d (red). The discovery limit for energy -|- time only readout is shown with the orange lines. The 
left (right) panel is for = 6 (100) GeV and detector set-up A (B), with low (high) mass and threshold. 


sional readout the decrease in sensitivity is larger and 
the plateau in the limit is reached for a larger head-tail 
energy threshold. Qualitatively similar behaviour occurs 
for the 100 GeV WIMP and the high mass and threshold 
detector. In this case the discovery limits flatten off to 
values 1.1 —1.2 below the energy only limit at a head-tail 
energy threshold of 60 keV. 

In Figure|^we show the discovery limits with and with¬ 
out sense recognition, as a function of WIMP mass. The 
factor by which the discovery limit changes without sense 
recognition is largest for light, < 0(20 GeV), WIMPs 
and a low threshold. The discovery limit achieved by 
a 3-d readout is still considerably lower than the non- 
directional limit however 1-d and 2-d readouts do suffer 
without sense recognition and are only marginally better 
than the non-directional limits, especially at high WIMP 
masses. 

In Figure [l0| we show (in similar fashion to Fig.|^ the 
evolution of the discovery limit now as a function of de¬ 
tector mass for = 6 and 100 GeV with and without 
sense discrimination. As in Fig. we see that the lack 
of sense recognition is most damaging in the 100 GeV 
WIMP case. This is particularly true for the 1-d and 2-d 
readouts where with no sense recognition there is only a 
factor of 1.1 and 1.2 improvement over a detector with 
no directional information at all and the evolution of the 
discovery limit suffers from the same saturation effect due 
to the similarity in the recoil distributions. In the 6 GeV 
WIMP case the discovery limits with no sense recognition 
continue to decrease past the saturation regime suffered 
by the non-directional limit. However there is still a re¬ 
duction in sensitivity by factors of 1.9, 2.8 and 8.9 for 
3-d, 2-d and 1-d readouts respectively compared to the 


limits with sense recognition. Interestingly, the discovery 
limit for 3-d readout with no sense recognition is slightly 
better than 1-d and 2-d readouts with sense recognition. 

Our main conclusion regarding sense recognition is 
that for discriminating between Solar neutrinos and low 
mass WIMPs, having 3-d readout with no method of de¬ 
termining sense is marginally preferable to 1-d or 2-d 
readout with sense determination. This is because the 
recoil distributions from low mass WIMPs and Solar neu¬ 
trinos are both anisotropic and have sufficient 3-d angular 
separation that they are still distinguishable even with¬ 
out recoil sense information. For the higher mass WIMPs 
this is not the case, and without sense recognition the ad¬ 
vantage of directionality is almost entirely lost, even in 
3-d. 


D. Angular resolution 

The final experimental limitation we study is that of 
finite angular resolution caused by the inaccuracy in the 
estimation of the “true” recoil direction. This is an in¬ 
herent difficulty faced by all directional detectors. For in¬ 
stance, directional detectors using low pressure gas TPCs 
suffer from straggling effects as the recoiling nucleus col¬ 
lides with other gas nuclei and more importantly from 
the diffusion of the primary electrons while drifting to¬ 
ward the anode m- Finite angular resolution will smear 
out the WIMP and Solar neutrino distributions in Fig.|^ 
making it more difficult to discriminate between the two. 
Since the minimum separation between the peak WIMP 
and neutrino directions is ^ 60°, an angular resolution 
better than this will likely be required to differentiate 
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FIG. 11: The discovery limit as a function of angular resolu¬ 
tion, (j-f for a detector with 3-d readout. The upper (lower) 
set of lines are for = 6 (100) GeV and detector set-up A 
(B), with low (high) mass and threshold. The dashed lines 
show the discovery limit using energy information only. 


FIG. 12: The discovery limit, as in Fig. as a function of 
WIMP mass with angular resolution, varying between 0° 
and 60° for a detector with 3-d readout. The upper (lower) set 
of lines are for detector set-up A (B), with low (high) mass 
and threshold. The dashed black lines show our discovery 
limit with energy information only and the black curve and 
shaded region shows the neutrino floor from Ref. m- 


between the WIMP and Solar neutrino distributions. 

Finite angular resolution results in a recoil in the di¬ 
rection r'(n(.) being reconstructed in the direction ^( 11 ^) 
with a probability distribution that has the form of a 
Gaussian smoothing kernel on a sphere [33 IHH] 


K{nr,n'^) 


1 

(27r)3/2cr..yerf(-\/2CT.y) ^ 



where 7 is the angle between the original and recon¬ 
structed directions. 


cos 7 = sin 0 sin 6' cos {(j) — (j)') + cos 9 cos 9 ', (25) 


in the co-ordinates dehned in Eq. |18| The measured 
directional recoil rate is then the convolution of the 
smoothing kernel with the original directional recoil rate 


d^i? 

dflrdEr 



d^i? 

dn(dEr 


{VL'^,Er)K{nr,n'^)dn'^. 


(26) 


The discovery limit as a function of angular resolution, 
is shown in Fig. 11 for = 6 and 100 GeV for our 
two example detector set-ups (Table [n| with 3-d readout. 
As expected, finite resolution makes it harder to discrimi¬ 
nate a 6 GeV WIMP from Solar neutrinos. The discovery 
limit is an order of magnitude weaker for = 30° than 
for perfect angular resolution, and for > 50° the limit 
is only marginally better than that obtained using en¬ 
ergy information only. For the heavier WIMP and the 
more massive detector, the discovery limit only has a 
slight change with increasing This is because the fi¬ 
nite angular resolution affects only the WIMP signal and 


not the isotropic background from atmospheric neutri¬ 
nos. However in this case the improvement afforded with 
directionality, even in the ideal case, is smaller. 

Figure shows the discovery limits for the two de¬ 
tector set-ups as a function of WIMP mass and angular 
resolution. Finite angular resolution significantly limits 
the ability of a low threshold directional detector to dis¬ 
criminate light, < 0(20 GeV), WIMPs from Solar 
neutrinos. The effects of hnite angular resolution are 
greatest for ^ 6 GeV, when the energy spectra of 
WIMPs and ®B neutrinos match one another. For the 
high threshold detector the reverse behaviour is observed. 
At higher WIMP masses the effect of increasing angular 
resolution is more apparent than at lower masses (< 12 
GeV), this is because the anisotropy of the recoil distri¬ 
bution decreases with increasing WIMP mass. 

The main conclusion of this sub-section is that angular 
resolution of order = 30° or better is required to ex¬ 
ploit the different directional signals of light WIMPs and 
Solar neutrinos. For angular resolutions larger than this 
there is little benefit from having the directional informa¬ 
tion at all as the Solar neutrino and WIMP signals are 
poorly resolved. For heavier WIMPs the neutrino floor 
can still be overcome even with angular resolutions up to 
60°. This is because the dipole asymmetry of the WIMP 
recoil distribution has a large dispersion and the effect 
of smearing due to finite angular resolution is less signif¬ 
icant. Therefore for light WIMPs probing cross-sections 
below the ®B neutrino floor requires good angular res¬ 
olution, however for the atmospheric neutrino floor the 
experimental limits can be competitive even with only 
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modest angular resolution. 

VI. CONCLUSIONS 

We have studied in detail how direct detection experi¬ 
ments with directional sensitivity can subtract the back¬ 
ground due to coherent neutrino-nucleus scattering and 
circumvent the so-called neutrino floor over a wide range 
of WIMP masses. In particular for light WIMPs (which 
have a similar recoil energy spectrum to ®B Solar neutri¬ 
nos) directionality would allow a ton-scale low threshold 
detector to be sensitive to cross-sections several orders of 
magnitude below the neutrino floor. We have also shown 
that experiments that can only measure 1 -d or 2 -d pro¬ 
jections of the recoil tracks can still discriminate WIMPs 
from neutrino backgrounds. 

Moving beyond ideal detectors, we studied the effects 
of finite angular resolution and limited sense recognition. 
The angular distributions of WIMP and Solar neutrino 
induced recoils are sufficiently different that for light 
WIMPs sense recognition is not crucial. The discovery 
limits are a factor of roughly two and ten worse without 
sense recognition for 3-d and 1-d readout respectively. 
However the discovery limits still improve strongly with 
increasing exposure. The discovery limit for 3-d read¬ 
out with no sense recognition is slightly better than 1 - 
d and 2-d readouts with sense recognition. For heavier 
WIMPs, however, sense recognition is required to dis¬ 
criminate WIMPs from the isotropic background from 
atmospheric neutrinos. Finally we found that if the angu¬ 
lar resolution is worse than of order thirty degrees, then 
it becomes signihcantly more difficult to discriminate be¬ 
tween light WIMPs and Solar neutrinos. Angular resolu¬ 
tion is less crucial for distinguishing heavier WIMPs from 
isotropic atmospheric neutrino (although in this case the 
improvement offered by an ideal directional detector is 
smaller). 

We have used Xenon as a target nucleus throughout, 
although we note that no directional dark matter de¬ 
tectors using this material currently exist (and the only 


proposed directional detection strategy using Xe has a 
1-d readout [55]). Using Xenon allows our results to be 
easily compared with previous studies of the neutrino 
floor HMm and simplifies the analysis as it is not nec¬ 
essary to consider the effects of multiple target nuclei, as 
is the case with CF 4 which is most commonly used in 
current low pressure gas TPCs jSlj- Moreover, multiple 
target experiments and their complementarity have al¬ 
ready been studied extensively in Ref. m- We have also 
not accounted for astrophysical uncertainties both in the 
values of parameters such as the escape speed and lab ve¬ 
locity but also in the shape of the velocity distribution. 
Ref. [T^ studied the effect on the light WIMP neutrino 
floor of a non-Maxwellian speed distribution and a distri¬ 
bution containing a stream. It found that the inclusion of 
time information allowed the additional uncertainty from 
the speed distribution to be overcome, and the neutrino 
floor suppressed. We expect that directionality would 
also help in a similar way. 

The results presented in this paper make a compelling 
case for the development of large directional dark 
matter detectors. If the results of the next generation 
of direct detection experiments lead the search to 
smaller WIMP-nucleon cross-sections, new techniques 
will need to be implemented to tackle the neutrino back¬ 
ground. We have shown that the use of directionality is 
a powerful way of doing this, even for non-ideal detectors. 
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